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We present a version of the Liischer-Weisz multilevel algorithm ideally suited for studying excited 
states of the QCD flux tube. While the original version achieved error reduction only in the 
temporal direction, the new algorithm reduces fluctuations in the sources as well. We report on 
the implementation of this algorithm as well as improvement over the older method. We also 
present first results, where we see a good agreement with theoretical predictions from bosonic 
string models. 
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1. Introduction 

Formation of a flux-tube between a quark and an antiquark in the QCD vacuum plays an 
important role in the description of the still unsolved phenomenon of quark confinement. At a large 
distance R between the quarks, it is believed that the QCD flux tube has quite a similar behavior to 
an oscillating bosonic string, largely independent of the internal gluonic degrees of freedom which 
make up the flux tube. 

This led to the formulation of effective string theories for the flux tube ^ and a lot of studies 
have dealt with the groundstate of the string [^. The simplest model describing the excited states of 
the string is that of the free bosonic string. A more realistic model is provided by the Nambu string 
theory. An important development by Arvis was the writing down of a qq potential assuming 
the flux tube to be described by a Nambu string. 

Nevertheless the Nambu theory still has the conformal anomaly in anything other than 26 
dimensions. This has led to several other proposals for effective string theories. LUscher and 
Weisz [Ql imposed an open- closed duality on the string partition function and obtained the Nambu 
string spectrum up to order Another effective string theory introduced by Polchinski and 

Strominger [||] consists of the most general terms at every order in 1 /R, which does not introduce 
the conformal anomaly. In this theory too the string spectrum agrees with Arvis up to corrections 
of order \/R^ |^. For a recent review see [0]. 

In this article we look at the excitation spectrum of the flux tube formed between a static quark 
and an antiquark in 3 dimensional SU(2) lattice gauge theory and compare our results with the 
predictions for open bosonic string spectra. For a similar study in the closed string sector see 

2. Sources and excited states 

In 2-1-1 dimensions, the energy states of the oscillating string can be classified by parity and 
charge-conjugation properties, (C,P), and one can identify the ground states in these channels with 
the lowest four energy states of the string: 

£'o^(+,+) E2^{-,-) E3^{-,+) (2.1) 

Our goal is now to measure these energy states on the lattice. 

The groundstate of the potential can be measured well with Polyakov loop correlators, which 
have the spectral representation {P{R,T) P{0,T)) = L~=o^n exp [—£"„(/?) T] . In this case the co- 
efficients bn are integers and En{R) are the energy states at quark separation R. Effective string 
prediction for ground state V{R) = Eo{R) is 

V{R) = -\im^ln{P{R,T)P{0,T))=Vo + aR-^{d-2)l + ... (2.2) 

where Vq is an unphysical constant and a is the string tension. The \/R term is the well known 
Luscher-term [^], reproduced by all effective string theories. Different models however give differ- 
ent predictions for the excitation spectrum and therefore one way to distinguish between different 
string-models is to measure the energy-differences between the excited states. 

The Polyakov loop correlators used to calculate the groundstate of the potential are bad esti- 
mators for the excited states. A better way of measuring the excited states is provided by Wilson 
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Figure 1: Left: Spatial lines that correspond to the sources on the lattice. Right: A Wilson loop with 
sources at the ends, that lie in the middle of the time sUces. The sUces with the solid lines are the time slices 
with fixed Unes during the sublattice updates. 
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Table 1: Left : Combination of the sources into (C,P) eigenstates. Right : Run parameters for the multilevel 
algorithm we used. 



loops. However normal Wilson loops, with straight spatial lines at the ends, will again project 
strongly on the ground state, but weakly on the excited states. In order to get a preferential cou- 
pling to the excited states, we use a set of wavefunctions at the ends of the loops, called sources. 
These sources correspond to spatial lines on the lattice, that replace the straight spatial lines of the 
Wilson loops. The set of sources we use is shown in figure ([l]). 

If one creates superpositions of these sources with well defined parity and charge conjugation. 



these 'channels' couple directly to the excited states according to equation [2.1| . The superpositions 
are shown in table (|l|). 



There have been other attempts to measure the excited states. Juge, Kuti and Momingstar [ |10| ] 
used asymmetric lattices with small physical temporal extent and a large number of basis states to 
measure the excited states in different field theories. In that respect our approach is complementary 
to theirs, as we use relatively few basis states, but use much larger physical temporal extents to 
exponentially suppress the contamination due to excited states. 

We can now calculate the expectation values of Wilson loops T) with channel K at the 

ends, which has the representation: 

W^(/?,r) = a^e-^'(«)^ (^l+/3^g-AiJ'^(«)7 + ...) (2.3) 
Thus we are able to measure the different energy states with the formula: 

+ corrections (2.4) 



E'^iR) = In 
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The ''corrections" in this expression are due to the higher states in the channel and, as can be seen 
from equation (23) are exponentially damped with T\ and T2. This is why one would like to go to 
Wilson loops with large time extends. 

The problem with large Wilson loops is of course the very small signal to noise ratio. 

Even for the numerical cheap SU (2) LGT in d=2+l, conventional methods do not work. One 
way of reducing the error is provided by the Liischer-Weisz algorithm [11], which leads to an 
exponential error reduction for the time transporters of the Wilson loops. Putting the sources on 
fixed time-slices of the lattice and using the sublattice updates to reduce the fluctuations of the 
time transporters produces good results and is practicable for loops with time extends up to about 
a Fermi, see [p^. 

However this is not the best way to use the algorithm because the fluctuations of the sources 
are not reduced with this method. 



3. The new method 

To achieve further error reduction we now move the sources from the fixed lines to the middle 
of the time sHces. Such a Wilson loop with sources at the ends is also shown in figure (|l|). 
Main advantages of such a procedure are: 

• The fluctuations of the sources are reduced by the sublattice updates as well. 

• One can use multihit on single links of the sources that leads to a further error reduction. 

It is also beneficial to use different numbers of sublattice updates for the time slices that contain 
the sources and the time slices that contain only time-transporters. In this way it is possible to 
choose parameters for the algorithm to optimize the noise to speed ratio for the single parts of the 
Wilson loops. 

Several tests show, that it is good to use more sublattice updates for the sources than for the 
time transporters for excited states. If one chooses the right parameters for the sources and the time 
slices, one is able to achieve an error reduction of ^(10) for the same computation time. 

A related algorithm was used by Kratochvila and de Forcrand to look at string breaking with 
Wilson loops 



4. First results 



In our first run we worked on a 24^-lattice for /? = 4 — 9 and on a 48^-lattice for /? = 10 — 12, 
with j8 = 5 (ro = 3.9536(3)), made 2000 total measurements using the scheme of sublattice updates 
shown in table (|l]). Compared to 1 12] we were able to increase the time extent of the loops from 
r = 2,4,6,8 tor = 4,6,8,12. 

The naive energies are calculated with the formula: 



Ea{R) 



1 



In 



Wa{R,T2) 
Wa(R,Ti) 



with T[ and T2=4 and 8 



(4.1) 



To get rid of the contamination due to the higher states in the channels, we use a fit to the form 



1 



In 



WaiR,T2) 
Wa{RJi) 



Ea{R) + 



1 

T2-n 



ae 



-bT, 



>(T2-Ti) 



(4.2) 
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Figure 2: Measured energy values and predictions of the Arvis potential: The grey crosses are the naive 
energy values and the blue points are the corrected ones. The red lines are predictions of the Arvis potential. 
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Table 2: Naive and corrected energies for the four lowest states in the different (C,P) channels. 

with the fit parameters E a, b,c to calculate the corrected energy £„. For the corrected energies E2 
and £3 we were only able to use Wilson loops with the time extends T = 4,6, 8. For the first two 
states we used Wilson loops with all time extents. The values we obtained are given in table(^. 

Even more interesting than the total energy values are the energy differences. Again to take 
the corrections into account, we use a three-parameter fit 



' in 



T2-n 



Wn{R,T2)Wo{R,T,) 
WniR,Ti)Wo{R,T2) 



-AEno + jr^^be-'^^ (l - ^'-''(^^-^')) , (4.3) 



with the parameters AEno,b,c. The corresponding values are shown in table (||). 

Finally in figure (^, we plot the total energies and in figure (^ the energy difference against 
the predictions of the Arvis potential [^], which gives the energy states 




En{R) = aRJl + —-. [n- — {d-2)\. (4.4) 
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Figure 3: The energy differences A£io, AE20 and AEj,q and the prediction of the Arvis potential: The 
grey crosses are again the naive differences and the blue points the measured values. The red lines are the 
predictions of the Arvis potential and the green lines are the predictions of the free bosonic string. 
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Table 3: Corrected energy differences of the excited states with the ground state. 

5. Conclusions 

We have discussed an algorithm that allows us to use the advantage of the Liischer-Weisz 
algorithm for the time transporters as well as for the sources and we presented our first results 
using the algorithm. We were able to go to much bigger Wilson loops than was possible previously. 
We see from figure (|3|) that at the same qq separation R the energy difference AEio is much closer 
to the predictions from Arvis than the energy differences AE20 and AE30, which agrees with the 
expactations in [Q]. Also the magnitude of the corrections seems to be larger for the higher states. 
Unfortunately the error reduction obtained with the naive states is not sufficient to get a signal for 
the third excited state beyond T = 8. At the moment it is not clear whether this is due to lack of 
statistics or poor choice of parameters in the multilevel algorithm. Investigation into this issue is 
ongoing. 
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